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ABSTRACT 


Title  of  Thesis:  Tropical  200  mb  Divergent  Circulations  and  Their 
Relation  to  Outgoing  Longwave  Radiation 

Frank  Sornatale,  Master  of  Science,  1988 

Thesis  directed  by:  Anandu  D.  Vernekar,  Professor, 

Department  of  Meteorology 

We  present  Annual  and  interannual  variability  of  the  200  mb  level 
mean  meridional  circulations  (Hadley  and  Ferrel  cells)  and  200  mb  tropical 

Cl*' i  '  ■  •v 

east-west  circulations  (Walker  circulations)  as  inferred  from  the  velocitv 

_ 

potential  field.  Also,  we  present  a  relationship^,  between  divergence  at  200 
mb  and  outgoing  longwave  radiation  (OLR). 

A  major  feature  of  the  200  mb  tropical  divergent  velocity  field  is  the 
dominance  of  the  western  Pacific  monsoon  region  as  the  source  for  meri¬ 
dional  and  zonal  overturning  motions.  Smaller,  but  significant  similar 
source  regions  are  found  over  South  and  Central  America  as  well  as  central 
Africa.  One  significant  variation  is  their  intra-seasonal  transitions  to  the 
summer  hemisphere.  The  northern  hemisphere  circulations  under-go  more 
variations  than  those  of  the  southern  hemisphere. 

->  Intra-annual  variations  of  the  equatorial  Hadley  cell  are  dominated  by 
the  variations  of  Asian  Monsoon  circulations.  The  seasonal  fluctuations  in 
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the  divergent  velocity  field  in  the  western  Pacific  contribute  to  a  strong 


annual  cycle  in  the  inter-annual  pattern  in  the  tropics. 


The  results  on  east-west  over-turning  motions  show  that  there  are 


three  cells  associated  with  the  three  ascending  branches  over  the  Asian  mon¬ 


soon,  South  and  Central  America,  and  Central  Africa.  The  most  dominant 


cell  is  the  Walker  cell  with  its  ascending  branch  over  the  western  Pacific 


(Asian  monsoon)  and  descending  branch  in  the  eastern  equatorial  Pacific. 


The  intensity  and  location  varies  within  the  year.  The  inter-annual  varia¬ 


bility  appears  to  be  associated  with  the  El  Nino  event. 


The  relationship  between  divergence  and  the  equivalent  blackbody  tem¬ 


perature  index,  derived  from  OLR  data,  is  significant  in  the  tropical  convec¬ 


tive  region.  The  relationship  is  weak  in  the  subsidence  region  away  from  the 


tropical  belt,  20  °  north  or  south. 
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Tropical  200  mb  Divergent  Circulations  and  Their  Relation  to 
Outgoing  Longwave  Radiation 

1.  Introduction 

Large-scale  tropical  irrotational  motions  develop  in  response  to 
differential  heating  in  the  atmosphere.  The  rotational  motions  are  generated 
by  the  earth’s  rotational  effects  through  the  divergence  term  (/+$•)  V-V* 
in  the  vorticity  equation  where  /  is  the  Coriolis  parameter  and  f  is  the 
relative  vorticity.  Because  /  is  proportional  to  the  sine  of  latitude,  the 
effects  of  this  term  increases  from  the  equator  to  the  poles.  While  large  scale 
motions  in  most  of  the  atmosphere  are  dominated  by  rotational  com¬ 
ponents,  the  irrotational  motions  are  important  in  the  tropics  where  the 
Coriolis  parameter  is  small. 

The  tropics,  which  cover  half  the  area  of  the  globe,  provides  a  large 
amount  of  the  energy  needed  to  drive  the  global  circulations.  The  energy,  in 
the  form  of  sensible  and  latent  heat,  is  transported  towards  the  equator  by 
the  trade  winds  in  the  lower  levels.  This  energy  is  then  lifted  and  con¬ 
densed,  thereby  providing  sensible  heat  and  potential  energy  which  is  tran¬ 
sported  poleward  by  the  circulations  of  the  upper  troposphere.  Deep 
cumulus  convection  is  considered  to  be  the  major  mechanism  in  the  vertical 
transport  of  low  level  energy  to  the  upper  troposphere. 

The  theories  of  axially  symmetric  (zonally  averaged)  components  of  the 
general  circulation  of  the  atmosphere  explain  the  maintenance  of  momen¬ 
tum,  heat  and  water  vapor  in  the  atmosphere.  Lorenz  (1967)  shows  the 
contributions  of  mean  meridional  (Hadley)  circulations  in  explaining  these 


maintenances.  The  mean  meridional  circulations,  by  definition,  exist  solely 
in  divergent  flow. 

Bjerknes  (1969)  and  Krishnamurti  (1971a, 1971b, 1973)  studied  the  axi¬ 
ally  asymmetric  components  of  tropical  circulations.  These  asymmetries  are 
the  result  of  variations  in  thermal  and  orographic  forcing  along  the  zonal 
plane.  One  component  of  these  circulations  is  the  zonal  over-turning 
motion  in  the  equatorial  Pacific  Ocean.  Bjerknes  referred  to  the  pressure 
and  temperature  oscillations  supporting  this  circulation  as  the  Southern 
Oscillation. 

The  Southern  Oscillation  (SO)  is  a  dipole  structure  in  the  pressure  pat¬ 
tern  with  centers  located  in  the  tropical  eastern  and  western  Pacific  Ocean. 
The  steepness  in  the  pressure  gradients  between  the  centers  is  a  measure  of 
intensity  of  these  over-tmning  circulations.  Bjerknes  (1969)  and  Arkin 
(1982)  have  shown  this  relationship  using  the  200  mb  tropical  wind  field 
and  the  Southern  Oscillation  Index  (SOI)  which  is  a  measure  of  surface  pres¬ 
sure  gradients.  They  also  show  an  association  between  the  SOI  and  the 
equatorial  Pacific  sea  surface  temperature  (SST). 

Bjerknes  (1969)  renewed  interest  in  this  zonal  circulation  when  he 
described  the  characteristics  of  this  thermally  driven  cell  he  called  the 
Walker  circulation.  Krishnamurti  (1973)  went  on  to  characterize  this  local 
Pacific  circulation  as  one  branch  of  a  more  extensive  east-west  circulation  he 
found  extending  across  the  entire  circumference  of  the  equatorial  belt.  The 
Walker  circulation  is  very  important  to  both  the  tropical  and  extra-tropical 
circulations.  It  is  maintained  by  rising  air  in  the  western  Pacific  'source’ 


region  of  warm  sea  surface  temperatures  and  descending  air  over  the  eastern 
Pacific  ’sink’  where  the  sea  surface  temperature  is  relatively  cold.  Warm  and 
moist  air  is  drawn  into  the  western  Pacific  low  level  convergence  center, 
rises  in  its  ascending  column,  and  empties  the  mass  of  air  into  an  upper- 
tropospheric  horizontal  divergence  center.  The  air  from  this  divergence 
center  becomes  part  of  the  200  mb  Hadley  and  Walker  circulations.  The 
Walker  branch  circulates  back  into  the  200  mb  convergence  center  support¬ 
ing  the  eastern  Pacific  ’sink’.  The  low-level  divergence  out  of  the  base  of 
this  ’sink’  increases  the  easterly  winds  which  in-turn  cools  the  SST  by 
upwelling  affects.  This  cooling  creates  a  stronger  SST  gradient  between  the 
source  and  sink  regions  which  is  associated  with  the  intensification  of  the 
Walker  circulation. 

In  mid-latitudes,  rotational  components  are  normally  an  order  of  mag¬ 
nitude  larger  than  irrotational  components.  In  the  tropics,  however,  the 
irrotational  component  is  approximately  of  the  same  order  as  rotational 
components.  The  tropical  upper-tropospheric  irrotational  circulation  drives 
both  the  global  meridional  and  the  tropical  zonal  circulations.  Several  stu¬ 
dies  of  this  circulation  include  those  by  Paegle  and  Paegle  (1978),  Cressman 
(1981),  and  Chang  and  Lau  (1982).  They  have  studied  the  variability  of 
tropical  circulations  with  particular  attention  to  sub-tropical  jet  stream 
enhancement  in  areas  of  upper-tropospheric  convergence.  In  a  recent  study, 
Lau  and  Boyle  (1987)  show  the  relationship  between  the  variability  of  the 
horizontal  divergence  field  and  the  strength  of  the  sub-tropica!  jet  stream. 

Studies  of  horizontal  divergence  have  been  limited  by  the  availability 
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of  accurate  observations.  Conventional  ground-based  observations  are  lim¬ 
ited  in  space  and  time  because  of  the  scarcity  of  observation  sites  over  oce¬ 
ans  and  in  the  third  world  countries  which  are  heavily  populate  the  tropical 
belt.  The  usefulness  and  accuracy  of  satellite-  derived  winds,  as  with  most 
indirect  measurements,  has  some  serious  problems  which  restrict  its  use. 
Hubert  and  Whitney  (1971)  identified  sever  possible  error  sources,  including: 
errors  due  to  uncertanties  of  cloud  heights,  errors  due  to  non-advective 
cloud  motions,  errors  of  measurement,  and  errors  in  tracking  cloud  targets. 
These  deficiencies  have  led  to  other  sources  for  wind  field  data. 

A  source  currently  under  investigation  is  the  use  of  data  from  NOAA 
polar  orbiting  satellites  which  provide  a  full  spatial  distribution  of  the  net 
outgoing  longwave  radiation  (OLR)  from  twice  daily  measurements.  The 
tops  of  deep  convective  clouds,  because  of  their  cold  cloud-top  temperatures 
and  high  emissivity,  emit  lower  radiation  than  the  earth’s  warm  surface 
under  clear  skies.  Therefore,  in  the  tropics,  OLR  fields  can  be  used  to  infer 
low-level  convection  supporting  deep  cumulus  cloud  structures.  Rupert  and 
Gray  (1976)  showed  a  large  positive  correlation  between  the  intensity  of 
convection  associated  with  tropical  upper-tropospheric  cloud  clusters  and 
horizontal  divergence.  From  these  observations  one  may  use  OLR  as  a 
proxy  for  low-level  convergence  or  upper-tropospheric  divergence. 

Since  1974  several  studies  have  tried  to  estimate  kinematic  fields  from 
OLR  data.  Julian  (1984)  and  Kasahara  et.  al.  (1987b)  have  utilized  OLR  to 
improve  the  accuracy  of  the  initial  field  used  in  numerical  weather  predic¬ 
tion  by  deriving  irrotational  components  of  wind  from  OLR  data.  These 


studies  were  attempting  to  establish  a  relationship  between  OLR  and  tropi¬ 
cal  divergent  winds  on  a  daily  basis.  Their  results  show  only  limited  success. 

The  purpose  of  this  study  is  to  understand  the  variability  of  the  tropi¬ 
cal  200  mb  divergent  circulation.  At  the  same  time  we  seek  a  relationship 
between  the  monthly  mean  OLR  and  200  mb  horizontal  divergent  winds. 
Towards  this  end,  we  analyze  the  intra-  and  inter-annual  variabilities  in  the 
200  mb  divergent  wind  field  and  relate  them  to  variabilities  in  the  OLR 
field.  Chapter  2  gives  a  brief  description  of  both  data  sets  and  their  origins. 
Variations  in  the  seasonally  averaged  data,  the  zonal  and  time-mean  zonal 
average  of  the  irrotational  winds’  meridional  component,  and  the  meridional 
and  time-mean  meridional  average  of  the  irrotational  winds’  zonal  com¬ 
ponent  are  shown  in  Chapters  3,  4,  and  5  respectively. 

In  Chapters  3  to  5  we  use  independent  OLR  data  in  the  form  of  OLR 
and  Equivalent  Black  Body  (EBB)  temperature  indexes,  as  a  proxy  for  con¬ 
vection  to  assess  our  results  from  the  velocity  potential  data.  Finally,  in 
Chapter  6  we  do  a  correlation  analysis  between  the  200  mb  divergence  and 
EBB  indexed  data  to  quantify  the  relationship  between  the  monthly  mean 


2.  Data 


a.  Velocity  Potentials 

The  velocity  potential  data  used  in  this  study  is  a  2.5 0  globally  grid- 
ded  field  provided  by  the  National  Meteorological  Center  (NMC),  Climate 
Analysis  Center  (CAC)  in  the  form  of  ninety-six  monthly  means  covering 
the  period  from  January  1979  to  December  1986.  Averaging  was  done  twice 
daily  for  the  00Z  and  12Z  runs  of  the  Global  Data  Assimilation  System 
(GDAS)  used  at  the  NMC. 

Prior  to  1978  the  GDAS  used  a  Hough  an?.lysis  technique  which  con¬ 
tained  no  divergent  motions.  On  22  September  1978  the  Hough  analysis 
was  replaced  by  an  Optimal  Interpolation  (01)  process  wl  ich  does  include 
divergence.  Kistler  and  Parish  (1978)  describe  the  three  major  components 
of  the  GDAS  in  detail.  They  consist  of  analysis,  initialization  and  predic¬ 
tion.  Prediction  is  only  indirectly  relevant  to  this  study  and  will  not  be  dis¬ 
cussed. 

Under  this  01  scheme,  the  objective  of  data  analysis  is  to  determine  the 
values  of  meteorological  variables  at  regularly  spaced  grid  points  from 
sparse  sets  of  randomly  spaced  observations.  Compounding  the  problem, 
each  observation  system  has  its  own  level  of  accuracy  in  reporting  the  data. 
Interpolation  methods  based  on  statistical  techniques  are  used  to  do  this 
interpolation  in  the  sigma  coordinate  frame  of  the  prediction  model.  Since 
meteorological  observations  can  not  provide  a  reasonably  populated  initial 
field  of  data,  the  6-hr  forecast  from  the  previous  numerical  model's  run  is 
used  as  the  first  guess  field  by  the  GDAS. 


Spatial  interpolation  are  used  to  match  randomly  located  observed 
data  to  the  same  globally  gridded  resolution  as  the  first  guess  field.  This 
guess  field  is  then  subtracted  from  the  observed  field  to  form  a  difference 
field  used  in  further  processing.  The  final  raw  analysis  is  produced  by 
linearly  combining  the  difference  and  the  first  guess  field.  The  combination 
is  weighted  by  quality  assessment  for  the  various  observations.  The  result 
field  is  a  new  (corrected)  analysis  field. 

Numerical  models  have  a  hard  time  utilizing  data  which  are  not  con¬ 
sistent  with  the  governing  equations  of  the  model.  Any  inconsistencies  with 
the  particular  model’s  mass  motion  balance  requirements  produces  inertial 
gravity  waves  which  the  model  might  treat  as  a  perturbation  in  the  flow. 
As  a  perturbation,  numerical  models  may  allow  the  very  small  gravity  wave 
to  grow  to  orders  of  magnitude  greater  than  its  actual  strength  in  the  real 
atmosphere.  Because  of  this  potential  to  distort  this  type  of  data,  it  is 
important  to  filter  out  these  inconsistencies  before  the  model  integration 
begins.  This  filtering  occurs  in  the  initialization  phase. 

The  initialization  process  consists  of  interpolation  techniques  and 
mass-motion  balance  requirements  which  have  undergone  several  changes 
between  January  1979  and  December  1986.  On  27  May  1980  the  GDAS 
converted  from  a  9-layer  global  gridpoint  domain  to  a  12-layer,  24-wave, 
global  spectra!  model.  This  change  reduced  the  problems  associated  with 
inconsistencies  in  the  data  described  above  and  according  to  Kistler  and 
Parish  (1982)  it  improved  wind  forecasts.  On  14  August  1982  an  isobaric 
multi-variate  01  (MOI)  scheme  was  implemented  which  expanded  the  model 


from  24  to  30  waves.  Dey  and  Morone  (1985)  give  the  full  details  of  the 
MOI  system  which  they  found  improved  the  upper-tropospheric  wind  fore¬ 
cast. 

Adiabatic  initialization  using  geostrophic  assumptions  were  used  until 
September  1984.  This  made  the  difference  field  consistent  with  the  GDAS 
first  guess  field  which  is  also  based  on  geostrophic  assumptions.  These 
assumptions  weaken  the  magnitude  of  divergence,  causing  a  dampening  of 
the  Hadley  and  Walker  circulations  produced  by  the  analysis.  In  addition, 
diabatic  heating  accumulates  during  the  first  24  hours  of  the  model’s 
integration.  The  lack  of  diabatic  heat  suppresses  vertical  motions  in  the 
barotropic  tropics  which  further  inhibits  the  formation  of  thermally  direct 
cells. 

In  September  1984  the  adiabatic  initialization  process  was  replaced  by 
a  diabatic  non-linear  normal  mode  initialization  (NNMI)  process.  Mohanty 
et.  al.  (1986)  reported  this  improves  the  vertical  consistency  of  the  model 
and  improves  the  accuracy  in  the  placement  of  tropical  convective  systems. 

In  May  of  1986  major  changes  to  the  GDAS  were  made  to  both  the 
analysis  procedures  and  to  the  global  forecast  model.  The  analysis  changes 
included  expanding  corrections  to  a  40-wave  rhomboidal  truncation,  increas¬ 
ing  the  number  of  observations  affecting  a  correction  from  20  to  33,  and 
including  TIROS  sounding  data  over  land  (above  100  mb).  At  the  same 
time  the  global  spectral  prediction  model  was  replaced  by  an  upgraded  spec¬ 
tral  model  used  to  make  medium  range  forecasts.  This  version  expands  the 
horizontal  resolution  from  30  to  40  waves  while  the  vertical  resolution  is 


increased  from  12  to  18  sigma  layers.  New  physical  parameterizations 
replaced  less  sophisticated  procedures  of  the  older  global  spectral  model. 
The  reader  is  directed  to  the  report  by  Dey  et.  al.  (1987)  for  the  specifics  of 
all  the  changes  implemented  at  this  time. 

The  procedure  for  separating  wind  vectors  into  rotational  and  irrota- 
tional  components  is  described  in  detail  by  Dey  and  Brown  (1976).  The 
'Helmholtz  theorem’  is  used  to  break  a  velocity  vector  into  its  components: 

=  (i) 

where  'I'  is  the  stream  function  (rotational),  x  1S  the  velocity  potential  (irro- 
tational),  and  V  is  the  horizontal  gradient  operator.  Dey  and  Brown  (1976) 
applied  this  theorem  with  finite  differencing  methods  to  compute  the  rota¬ 
tional  and  irrotational  wind  components: 

f  -  V2  ¥  (2) 

D  =  V2  X  (3) 

where  D  is  the  divergence  and  f  is  the  relative  vorticity.  In  this  study  we 
have  chosen  the  convention 

-D  -  V2  X  (4) 

This  convention  reverses  the  signs  of  the  original  x  data  from  NMC,  which 
followed  the  conventions  of  Dey  and  Brown.  In  our  study  positive  values 
represent  divergence  and  negatives  represent  convergence.  The  irrotational 
wind  component  was  computed  from  the  velocity  potential  field  using  the 


following  finite  differences: 
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We  compared  our  calculations  with  those  of  Dey  and  Brown  and  found  they 
were  identical. 
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b.  Outgoing  Longwave  Radiation  (OLR) 

The  OLR  data  used  in  this  study  are  total  flux  estimates  obtained 
from  the  NOAA  polar-orbiting  satellite  series  which  are  low  orbiting,  sun- 
synchronous  satellites,  operating  since  1974.  The  data  were  digitized  on  a 
2.5  °  global  grid  and  consist  of  mean  monthly  values  calculated  from  twice- 
daily  observations  collected  from  January  1979  to  December  1986. 

The  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  aboard  the 
NOAA  Polar  Orbiters  measures  the  earth’s  terrestrial  radiation  budget  in 
the  10  -  12.5  /z m  window.  These  narrow  band  measurements  are  used  to 
estimate  the  total  outgoing  flux  over  the  full  spectrum  by  applying  empiri¬ 
cal  relationships  developed  by  Abel  and  Gruber  (1979).  The  collection  and 
interpolation  process  has  several  areas  of  bias  which  we  will  summarize  here. 

One  bias  exhibited  by  the  satellites  is  discussed  by  Gruber  and  Kruger 
(1984)  and  concerns  changes  in  the  equatorial  crossing  times  of  different 
satellites.  During  this  period  crossing  times  varied  from  the  early  0230-1430 
LST  crossings  of  NOAA7  to  the  late  0730-1930  LST  crossings  of  NOAA-6. 


These  variations  would  be  especially  noticeable  over  the  sub-tropical  desert 
regions  where  a  crossing  during  the  peak  heating  of  afternoon  (1430  LST) 
would  cause  a  positive  bias  in  the  flux  measurement.  Because  this  study 
uses  monthly  means  from  twice  daily  averaged  measurements,  we  assume 
the  affects  of  these  variations  should  not  be  significant. 

Another  bias  is  in  the  over-estimate  of  flux  values.  One  reason  is  dis¬ 
cussed  by  Ramanathan  and  Briegleb  (1980)  based  on  the  effects  of  cirrus 
clouds.  Cirrus  clouds  are  semi-transparent  in  the  infra-red  region  which 
lowers  their  emissivity  which  causes  an  over-estimate  in  the  flux  measure¬ 
ments.  Related  over-estimates  investigated  by  Ellingson  and  Ferraro  (1983) 
and  Ohring  et.  al.  (1984)  amounted  to  a  10  W  m~2  from  middle  and  high 
partial  cloudiness.  In  this  study  we  use  monthly  mean  values  which  we 
assume  reduces  the  occurrence  of  false  readings  from  these  over-estimates. 

In  this  study  we  use  OLR  indexes  to  infer  the  location  of  low-level  con¬ 
vergence  and  their  associated  upper-tropospheric  divergence  centers.  The 
procedures  we  use  provide  an  index  from  which  we  can  infer  the  intensity 
and  locations  of  convection  and  the  kinematic  motions  associated  with  it. 
Similar  inferences  between  maxima  in  the  OLR  data  and  upper-tropospheric 
convergence  centers  are  not  as  reasonable  (see  Kasahara  et.  al.,  1987b). 
Clear  skies  implied  by  the  maxima  can  be  caused  by  several  mechanisms 
including:  dry  convection  found  over  deserts,  upper-tropospheric  subsidence, 
or  neutral  activity  (neither  convection  or  subsidence). 

Furthermore,  the  link  between  convection  and  OLR  is  only  valid  in  the 
tropics  where  vertical  motions  are  proportional  to  the  diabatic  heating.  Fol- 


lowing  the  scaling  arguments  of  Holton  (1979),  both  the  local  temperature 
variation  and  the  temperature  advection  term  in  the  tropics  are  negligible, 

so  we  can  reduce  the  thermodynamic  equation: 

' 

4-+V-V  T-oSp  =  -3-  (7) 

J  CP 

to 


“  -  -rr  (8) 

*p  cp 

where  Sp  is  the  static  stability  term,  q  the  rate  of  diabatic  heating,  and  u> 
the  vertical-p  velocity.  In  tropical  precipitating  systems  the  vertical  velocity 
is  3  cm  s-1,  which  is  an  order  of  magnitude  larger  than  outside  of  the 
precipitation  regions.  In  the  extra-tropics  the  temperature  advection  can  not 
be  neglected  hence  the  relationship  we  have  established  here  between  OLR 
and  kinematic  motion  is  not  valid.  Therefore,  our  discussions  are  generally 
confined  to  the  tropical  belt  between  30 0  N  and  30  °  S. 


We  use  an  indexing  method  to  highlight  convection  inferred  by  the 
lowest  flux  measurements  and  the  clear  skies  by  the  highest  fluxes.  By  sub¬ 
tracting  OLR  flux  from  the  median  tropical  flux  (  240  W/m2)  we  can 
highlight  maxima  and  minima  in  the  OLR  field. 


OLR  Index  (OLRI)  -  240  W/m2  -  OLR  Flux  (9) 

In  this  index,  sign  conventions  have  convection  inferred  in  the  maxima  and 

clear  skies  in  the  minima.  This  conforms  to  conventions  for  200  mb  diver¬ 
gence  derived  from  velocity  potential  data. 


c.  Equivalent  Black  Body  (EBB)  Temperature  (Te) 

Several  studies,  including  those  of  Julian  (1984)  and  Kasahara  et.  al. 
(1987b),  use  Equivalent  Blackbody  (EBB)  temperature  (Te)  to  study  the 
relationship  between  divergence  and  OLR.  The  EBB  temperature  is  com¬ 
puted  using  the  ’Stephan-Boltzman  law’: 


where  <r  is  the  Stephan-Boltzman  constant  (5.6693  X  10  8  W  m  2  K  4). 
We  compute  the  EBBI  by  the  following  method: 

EBB  temperature  index  (EBBI)  =  258°K  —  Te  (11) 

The  258  °K  threshold  was  found  by  Kasahara  (1987a)  to  be  the  Te  separat¬ 
ing  rising  and  sinking  motions.  This  method  conforms  the  EBBI  sign  con¬ 
vention  to  those  described  for  OLRI  above. 


3.  Mean  Seasonal  Velocity  Potential  and  OLR 

In  this  chapter  we  present  the  mean  seasonal  variations  of  velocity 
potentials  and  Outgoing  Longwave  Radiation  (OLR)  for  the  8  years  from 
January  1979  to  December  1986.  As  described  in  Chapter  2,  the  velocity 
potential  field  was  affected  by  four  changes  in  analysis  procedures.  In 
proceeding  with  our  computations  we  assumed  these  changes  only  affect  the 
amplitudes  of  the  data  but  not  the  phases. 

In  this  chapter  maxima  in  the  velocity  potential,  OLRI,  and  EBBI 
fields,  denote  low  level  convergence  and  upper-tropospheric  horizontal  diver¬ 
gence.  Irrotational  velocity  vectors  computed  from  the  velocity  potential 
field  (see  Appendix)  have  been  superposed  on  figures  in  this  section.  The 
seasons  we  refer  to  are  based  on  the  northern  hemisphere  as  follows:  winter 
from  December  to  February  (DJF),  spring  from  March  to  May  (MAM), 
summer  from  June  to  August  (JJA),  and  fall  from  September  to  November 
(SON). 

Figure  1  shows  the  mean  DJF  distribution  for  the  velocity  potential 
and  OLRI  fields.  The  major  feature  in  velocity  potential  is  the  High-Low 
dipole  found  in  the  tropical  belt  which  dominates  the  global  pattern.  The 
irrotational  vectors  indicate  the  Borneo-Indonesia  maritime  continent  in  the 
southwest  Pacific,  which  we  will  refer  to  as  the  Pacific  Maritime  Continent 
(PMC),  is  the  central  location  for  the  major  north-south  (Hadley)  and  east- 
west  overturning  motions.  The  activity  in  the  velocity  potential  field  com¬ 
pares  favorably  with  the  OLRI  field  which  shows  the  major  convection 
region  in  the  vicinity  of  the  PMC.  Smaller  convection  areas  are  also  found 


over  south-central  Africa  and  central  South  America. 

Figure  2  magnifies  the  tropical  features  of  the  EBBI  field  during  DJF 
with  the  irrotational  vectors  superposed.  The  dominant  center  of  inferred 
convection  from  the  EBBI  field  is  over  the  PMC  with  a  narrow  zone  extend¬ 
ing  to  the  east  along  the  northwest  coast  of  Australia  and  then  to  the 
southeast.  The  strongest  velocity  vectors  originate  from  the  PMC  and  their 
magnitude  indicates  the  strength  of  the  north-south  overturning  motions. 
The  strongest  200  mb  divergent  winds  flow  towards  the  northwest  where 
Figure  1  shows  convergence  into  the  200  mb  sink  over  Tibet.  This  conver¬ 
gence  should  be  accompanied  by  a  strong  local  minima  in  the  OLRI  field. 
Instead,  a  local  maxima  appears  because  of  the  high  altitude.  The  surface 
temperatures  in  the  Himalayas  under  clear  skies  are  very  cold  which  cause 
low  OLR  flux  values.  This  is  one  example  where  topographic  effects  negate 
our  general  assumption  inferring  convection  from  minima  in  OLR  flux. 

Another  major  feature  is  the  east-west  circulations  extending  from  the 
PMC.  One  branch  flows  towards  the  upper-tropospheric  sink  in  the  south 
Pacific  convergence  Zone  (SPCZ)  and  another  to  the  sink  off  the  east  coast 
of  Africa.  A  second  east-west  flow  extends  from  the  South  American  con¬ 
vective  center  to  the  west  coast  of  Africa.  These  east-west  overturning 
motions  on  the  zonal  plane  will  be  looked  at  more  closely  in  section  5. 

Several  interesting  features  are  observed  in  Figure  2  around  Central 
and  South  America.  The  EBBI  field  infers  convection  over  central  South 
America  which  is  also  reflected  in  the  irrotational  velocity  vectors.  The  irro¬ 
tational  velocity  vectors  also  show  strong  convergence  zones  extending  from 
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the  Caribbean  Sea  to  the  western  Atlantic  Ocean,  and  off  the  coast  of 
Venezuela.  The  EBBI  field  infers  convection  over  south-central  Africa  and 
clear  skies  over  north-central  Africa.  The  irrotational  velocity  vectors  show 
a  convergence  sink  at  200  mb  off  the  southeast  African  coast. 

We  compared  our  results  with  a  similar  study  by  Lau  and  Boyle  (1987) 
which  showed  the  mean  velocity  potential  and  irrotational  wind  vectors  for 
the  9  DJFs  from  1974/75  to  1982/83.  The  general  patterns  in  the  two  stu¬ 
dies  agreed,  especially  the  divergence  pattern  associated  with  the  PMC. 
There  are  some  differences  however,  mainly  in  the  intensity  of  the  irrota¬ 
tional  wind  vectors  over  South  America,  the  Caribbean,  and  the  western 
Pacific.  There  was  also  some  disagreement  on  the  spatial  placement  of 
several  minor  centers  of  minima  and  maxima  values.  We  attribute  these 
disagreements  to  differences  in  the  data  sets  between  the  two  studies. 

Figures  3  and  4  compare  the  data  fields  for  MAM.  The  maximum  velo¬ 
city  potential  in  the  PMC  decrease  from  9.2  m/s2  in  DJF  to  7.25  m/s2  in 
MAM.  Accompanying  the  decrease  are  weaker  irrotational  winds  which 
decrease  from  4.4  m/s  to  2.37  m/s.  The  irrotational  winds  are  more  sym¬ 
metric  about  the  equator  than  they  were  during  DJF  and  the  dominant  flow 
pattern  is  now  directed  to  the  south  and  southeast  out  of  the  PMC.  The 
core  of  the  200  mb  divergence  shifted  to  the  north  and  a  minima  in  velocity 
potential  appears  off  the  west  coast  of  South  America. 

The  OLRI  pattern  shows  a  northward  drift  of  the  PMC  convection 
while  a  decrease  in  velocity  potential  gradients  is  apparent.  The  convection 
implied  by  OLRI  maxima  over  South  America  and  Africa  covers  less  area 


The  details  of  Figure  4  show  a  realignment  and  intensity  drop-off  in 
convection  implied  by  the  decrease  in  the  PMC  EBBI  from  18 4  K  in  DJF  to 
15 4  K  in  MAM.  The  lower  index  value  shows  the  top  of  the  convection  is  at 
a  lower  level  in  the  atmosphere  and  radiates  at  a  warmer  temperature  than 
those  with  higher  index  values.  From  the  weaker  convection  we  infer  that 
supporting  kinematic  motions  are  also  weaker,  specifically  the  horizontal 
divergence  at  the  outflow  of  the  convection.  The  200  mb  divergence, 
reflected  by  the  irrotational  wind  field,  has  moved  to  the  north  and  the 
magnitude  of  the  flow  has  decreased  since  DJF.  The  local  centers  over 
South  America  and  Africa  have  also  undergone  some  realignment  since  DJF 
having  moved  approximately  5  0  to  the  north.  There  is  stronger  convergence 
off  the  west  coast  of  South  America  than  in  DJF  but  the  convergence  center 
off  the  northwestern  coast  has  decreased  considerably.  The  EBBI  field 
implies  strong  convection  in  central  Africa  and  convergence  in  northeast 
Africa.  These  patterns  are  not  clearly  visible  in  the  irrotational  wind  field. 

The  JJA  patterns  are  shown  in  figures  5  and  6.  The  200  mb  divergence 
center  over  the  PMC  continues  to  move  northward  to  approximately  10  0  N. 
The  OLRI  inferred  convection  supporting  this  circulation  appears  to  be 
bordering  the  equator  but  extending  away  from  the  PMC  toward  the 
Southeast  Asian  peninsula.  The  maximum  horizontal  divergent  velocity  vec¬ 
tor  is  4.98  m/s  which  is  the  stronger  than  during  any  other  season. 

During  the  transition  from  MAM  to  JJA  the  effects  of  the  summer 
monsoon  become  evident.  Over  the  latter  part  of  the  MAM  season,  when 


south  Asia  is  warming  and  the  Southern  Hemisphere  is  cooling,  the  pressure 
gradients  force  air  northwards  forming  a  Low  pressure  system  over 
southeast  Asia.  This  is  accompanied  by  increased  rainfall  over  Burma  and 
rising  temperatures  over  India.  The  cycle  breaks  in  late  May  to  early  June, 
with  intense  storm  activity  moving  to  the  northeast  from  Burma.  This 
monsoon  cycle  is  visible  in  our  figures  even  though  the  data  is  a  three- 
month  average. 

Figure  5  shows  that  the  summer  monsoon  has  moved  northwards 
across  the  equator  and  broadened  considerably  since  MAM.  The  convection 
now  extends  from  the  PMC  northeastward  encompassing  the  Indian  Ocean, 
Southeast  Asia,  and  the  Indian  subcontinent.  The  narrow  zone  of  convec¬ 
tion  in  the  central  Pacific  is  seen  in  the  EBBI  pattern  paralleling  the  equa¬ 
tor  at  10 0  N  latitude,  which  is  the  preferred  location  of  the  inter-tropical 
convergence  zone  (ITCZ)  during  JJA.  The  EBBI  pattern  at  30  °  N  150  0  W 
indicates  an  area  of  clear  skies  which  is  j  sistent  through  all  seasons  but 
has  its  most  extreme  values  in  this  JJA  pattern. 

Two  east-west  overturning  circulations  come  from  the  summer  mon¬ 
soon  convection.  One  branch  is  in  the  NH  moving  toward  the  west  coast  of 
Mexico,  and  a  second  is  moving  in  the  SH  toward  the  west  coast  of  South 
America.  Very  strong  zonal  irrotational  winds  flow  from  the  summer  mon¬ 
soon  towards  the  west  coast  of  Africa. 

The  northwestward  shift  in  convection  is  very  noticeable  in  the  EBBI 
patterns  over  South  America.  The  inferred  convection  has  moved  from  cen¬ 
tral  South  America  in  DJF  to  the  Central  American  isthmus  in  JJA.  In 


addition,  it  has  developed  an  elongated  shape  from  Venezuela  to  the  coastal 
waters  of  Baja  Mexico  in  alignment  with  this  land  mass.  The  divergent 
velocity  vectors  show  support  for  this  inferred  convection  and  the  clear  sky 
region  near  the  SPCZ. 

In  MAM  the  elongated  maxima  in  the  EBBI  pattern  over  Africa  has 
split  into  two  cells.  Divergent  velocity  vectors  flow  from  this  convection  to 
the  upper-tropospheric  convergent  zone  off  the  southwestern  coast.  There  is 
also  support  in  the  irrotational  winds  for  the  central  African  convection  and 
the  200  mb  convergence  zone  in  northeast  Africa  and  Pakistan. 

The  final  figures  (7  and  8)  show  the  mean  SON  season.  There  is  a  gra¬ 
dual  shift  in  the  activity  back  to  the  south  and  a  weakening  in  magnitudes 
of  both  fields  from  the  JJA  maxima.  Figure  8  clearly  shows  a  long  streak  of 
inferred  convection  across  the  Pacific  in  the  position  of  the  Inter-Tropical 
Convergence  Zone  (ITCZ).  The  clear  areas  in  the  central  and  eastern  Pacific 
are  weaker  but  cover  a  very  large  area  parallelling  the  ITCZ.  The  conver¬ 
gence  off  the  West  coast  of  South  America  is  stronger  than  any  other  sea¬ 
son.  The  inferred  convection  in  the  EBBI  field  is  less  intense  than  JJA  but 
more  centralized  over  Venezuela  feeding  the  upper-tropospheric  convergence 


zone  off  the  coast. 

The  onset  of  the  cold  season  over  Siberia  forces  the  surface  winds  we 
discussed  in  JJA  to  reverse  direction  from  ’southerlies’  to  ’northerlies'.  This 
increases  the  convergence  into  the  PMC  and  signals  the  start  of  the  winter 
monsoon  season.  Over  Africa  the  EBBI  pattern  shows  the  western  convec¬ 


tive  cell  moves  off-shore  to  the  central  Atlantic.  The  divergent  winds  sup- 
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port  the  inferred  convection  and  subsidence  in  the  EBBI  field  over  central 
and  northern  Africa. 

We  compared  the  patterns  of  Figures  1-8  with  the  January,  April, 
July,  and  October  figures  in  Arkin  et.  al.  (1986)  which  was  a  similar  study 
of  monthly  mean  velocity  potentials  from  October  1978  to  September  1983. 
The  velocity  potential  distributions  are  identical  except  for  sign.  We  use  the 
convention  that  irrotational  wind  blows  from  High  velocity  potentials  to 
low  velocity  potentials,  Arkin  et.  al.  use  the  opposite  convention. 

Several  characteristics  of  the  intra-annual  cycle  appear  in  these  figures. 
The  first  is  the  dominance  of  the  PMC  divergence  center  in  governing  the 
global  Hadley  regime.  Second  is  the  seasonal  variation  in  the  location  and 
intensity  of  this  center  in  alignment  with  the  winter  and  summer  monsoons. 
It  is  evident  from  the  horizontal  irrotational  wind  field  that  this  western 
Pacific  High  is  the  source  for  the  tropical  Hadley  branch  of  the  general  cir¬ 
culation.  As  we  shall  see  in  the  next  section,  the  seasonal  variation  of  the 
Hadley  cell  is  dominated  by  the  seasonal  north-south  movement  of  this  con¬ 
vection. 

The  200  mb  PMC  outflow  is  stronger  to  the  west  than  to  the  east.  We 
also  observed  differences  in  hemispheric  patterns.  Some  of  these  differences 
are  related  to  the  zonaliy  asymmetric  distribution  of  land/sea  mass  about 
the  equator.  The  southern  hemisphere  (SH)  patterns  are  more  persistent 
because  of  the  stable  heating  gradients  associated  with  its  large  water  mass. 
In  the  northern  hemisphere,  large  land  masses  contribute  to  heating 


differentials  between  the  land/sea  topography.  This  causes  substantial 


fluctuations  in  the  circulations. 

The  patterns  in  the  southern  hemisphere  show  only  the  slightest  varia¬ 
tions  but  magnitudes  do  vary  between  seasons.  The  western  flow  from  the 
PMC  source  converges  into  the  upper-tropospheric  sink  off  the  east  African 
coast  and  the  eastern  flow  converges  into  the  SPCZ  sink.  These  ’westerlies’ 
reorganize  off  the  east  coast  of  South  America  and  converge  into  the  sink  off 
the  southwest  coast  of  Africa. 

The  EBBI  pattern  infers  a  continuous  convective  belt  see-sawing  along 
the  equator  which  is  associated  with  the  seasonal  changes.  Clear  sky  returns 
bound  this  convective  band  to  the  north  and  south.  The  dominant  convec¬ 
tion  is  in  the  vicinity  of  the  PMC  but  varies  according  to  the  monsoon 
activity  in  this  area.  The  most  intense  activity  is  during  the  summer  mon¬ 
soon  when  convection  spreads  from  the  PMC  to  the  Indian  subcontinent. 
This  variation  in  intensity  parallels  changes  in  the  intensity  of  the  200  mb 
irrotational  wind  field.  From  a  qualitative  view,  it  appears  velocity  poten¬ 
tial  and  OLR  data  agree  on  the  inferred  spatial  placement  of  upper- 
tropospheric  tropical  divergence  and  convergence  centers. 
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Figure  1.  DJF,  (a)  Mean  200  mb  velocity  potential  and  irrotational  velocity 
vectors,,  contour  interval  2  x  106  m2  s~2.  (b)  OLRI,  contour  interval  10  w 
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Figure  5.  As  in  Figure  1  except  for  JJA. 
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4.  Mean  Meridional  Circulation 

In  Chapter  3  we  showed  the  seasonal  variations  of  the  velocity  poten¬ 
tial  field  and  the  200  mb  divergent  circulation  associated  with  overturning 
motions  in  the  meridional  and  zonal  planes.  In  this  section  we  look  at  the 
intra-  and  inter-annual  variations  of  zonal  mean  meridional  velocity,  200 
mb  divergence,  and  zonal  mean  EBBI. 

Figure  9a  shows  the  intra-annual  variations  of  zonal  mean  meridional 
velocity.  The  positive  (negative)  values  represent  northerly  (southerly)  flow. 
The  arrows  which  have  been  superposed  on  some  of  these  figures  are  meant 
to  highlight  the  inferred  direction  of  the  flow  and  are  not  related  to  their 
intensities. 

The  zero  contour  near  the  equator  in  Figure  9a  is  the  200  mb  diver¬ 
gence  line  which  represents  the  ascending  branch  of  the  equatorial  Hadley 
cell.  During  DJF  this  branch  is  found  at  20  “  S  and  the  inter-hemisphere 
flow  is  directed  from  the  southern  hemisphere  (SH)  to  the  northern  hemi¬ 
sphere  (NH).  This  pattern  reverses  by  JJA  when  the  ascending  branch  is 
found  at  30  °N  and  the  inter-hemisphere  flow  is  now  from  NH  to  SH.  The 
descending  branches  of  the  equatorial  Hadley  cell  are  near  30 0  N  and  30  °  S 
during  DJF.  In  the  SH  it  remains  near  30  0  S  throughout  the  year  but  in  the 
NH  it  is  in  the  vicinity  of  50 0  N  in  JJA. 

Oort  and  Rasmusson  (1970)  studied  the  inter-annual  variations  of 
mean  meridional  circulations  based  on  the  5-year  period  from  May  1958  to 
April  1963.  Our  results  generally  agree  with  theirs  but  there  are  some 
significant  differences.  The  major  difference  is  that  our  magnitudes  were 


considerably  smaller  than  theirs  particularly  in  the  northern  hemisphere 
DJF  where  their  values  were  three  times  larger  than  ours.  They  show  the 
maxima  in  the  JJA  ascending  branch  in  August  while  our  maxima  occurs  in 
July  and  is  found  10  *  further  south.  Another  difference  is  in  the  NH  above 
30  *  N  where  the  lip  of  northerly  flow  extending  to  35  0  N  during  JJA  is  not 
found  in  their  results. 

Figure  9b  is  the  inter-annual  cycle  of  the  zonal  mean  meridional  velo¬ 
city.  The  vertical  lines  superposed  on  this,  and  Figure  11a,  show  the  month 
when  analysis  changes  to  the  GDAS  were  implemented.  These  partitions 
highlight  the  major  weakness  of  this  data  set  which  is  the  inhomogeneity 
over  the  8  year  period  caused  by  the  changes.  At  the  same  time,  these 
figures  support  the  assumption  stated  in  Section  3  that  phases  remain  con¬ 
stant  through  the  data  changes. 

The  figure  shows  a  strong  cellular  cycle  in  the  tropical  belt  from  30  0  N 
to  30  °  S  where  northerly  winds  are  replaced  by  southerly  winds  with  6 
month  period.  This  periodicity  is  associated  with  the  seasonal  fluctuations 
of  the  equatorial  Hadley  cell.  This  cycle  remains  both  spatially  and  tem¬ 
porally  consistent  over  the  8  years  with  maximum  southerly  winds  at  10  0  N 
during  DJF  and  maximum  northerly  winds  at  5  °  S  during  JJA.  The  only 
variations  are  in  the  magnitudes,  particularly  in  1986  after  the  last  GDAS 
changes  were  implemented.  In  JJA  of  1986  the  northerly  wind  has  a  magni¬ 
tude  of  -4.25  m/s,  compared  to  values  in  the  range  of  1.9  m/s  to  1.6  m/s 
from  1980  to  1984.  In  spite  of  magnitude  variations,  the  annual  cycle 
remains  very  consistent  over  the  entire  8  years. 


The  oscillations  north  of  30  °N  are  the  variations  associated  with  the 
NH  Ferrel  cell.  This  cell  shows  similar  oscillations  to  those  of  the  equatorial 
Hadley  cell  but  the  SH  Ferrel  observed  south  of  30  °S  shows  hardly  any 
variations.  The  one  exception  is  during  DJF  of  1983  and  1984  when  the 
northward  flow  pattern  extends  across  20  °  S  bridging  through  to  the  des¬ 
cending  motions  of  the  SH  Ferrel  cell  at  30  °  S.  In  the  divergence,  as  we  will 
see,  there  is  also  some  discontinuity  during  these  periods.  These  variations 
are  possibly  related  to  the  1982-83  El  Nino  event. 

Figure  10a  shows  the  intra-annual  cycle  of  the  mean  zonal  divergence 
and  10b  the  mean  zonal  EBBI  field.  We  have  marked  the  maxima  in  to 
show  lines  of  divergence  and  convergence.  Maxima  in  equatorial  divergence 
begin  at  10  °S  '...  January  and  moves  to  10  °N  by  May.  It  remains  here 
until  November  before  moving  back  to  10  *  S.  The  SH  convergence  line 
varies  slightly  from  35  °  S  in  January  to  25  °  S  in  July.  The  NH  convergence 
maxima  line  appears  near  30  °N  in  January  and  remains  there  through 
June. 

The  mean  zonal  EBBI  field  (10b)  generally  agrees  with  the  divergence 
field.  The  line  of  inferred  convection  at  10°  S  in  January  moves  to  5  0  N  by 
May  and  then  back  to  10  °  S  by  December.  The  clear  skies  in  the  NH  pat¬ 
tern  is  represented  along  a  line  starting  at  10  °N  in  January,  moves  to 
20°  N  by  May,  then  back  to  10°  N  by  December.  This  pattern  of  equatorial 
convection,  undergoing  seasonal  cross-equatorial  transitions,  is  found  in  all 
of  the  figures. 

Figure  11  shows  a  comparison  between  the  mean  zonal  divergence  and 


the  EBBI  field.  The  divergence  has  a  High-Low  cellular  structure  across  the 
equator  throughout  the  8  years.  This  pattern  is  in  phase  with  annual  sea¬ 
sonal  variations.  Large  magnitude  are  again  observed  after  the  last  GDAS 
change  in  1986  was  implemented.  In  this  figure  a  maximum  infers  upper- 
tropospheric  horizontal  divergence. 

The  clearest  association  between  these  fields  is  the  wavy  pattern  along 
the  equator  between  20  °  N  to  20  °  S  over  the  full  8  year  cycle.  The  EBBI 
pattern  is  more  regular  in  appearance  and  maintains  its’  shape  over  the  8 
years.  The  divergence  field  has  a  slightly  jagged  appearance  with  slight 
breaks  in  the  SH  pattern  during  DJF  of  1983  and  1984.  This  was  the  same 
type  of  variation  we  observed  in  the  meridional  pattern  (9b). 

These  figures  infer  bands  of  convergence  (lla)  and  clear  skies  (lib) 
between  15  °  N  and  40  °  N  in  the  NH  and  5  °  S  and  30  °  S  in  the  SH.  This 
implied  subsidence  region  borders  the  divergence  wave.  The  only  disparity  is 
in  the  NH  of  lla  where  branches  of  divergence  extend  across  the  NH  band 
during  JJA. 

In  this  section  we  observed  the  intra-  and  inter-annual  variability  in 
the  zonal  mean  meridional  velocity,  zonal  mean  divergence,  and  mean  EBBI 
field.  The  same  general  patterns  were  supported  in  the  divergence  and  the 
EBBI  field  which  implies  that  preferred  spatial  and  temporal  zones  of  200 
mb  divergence  are  supported  by  two  independent  data  sets.  The  seasonal 
meridional  transition  that  was  so  prevalent  in  Chapter  3,  was  supported  by 
both  intra-  and  inter-annual  fields.  These  patterns  support  the  global  struc¬ 
ture  of  the  tropica!  Hadley  circulation  and  its  characteristic  of  cross- 


equatorial  flow  patterns.  Oort  and  Rasmusson  (1970)  showed  similar  pat¬ 
terns  but  some  discrepancies,  especially  in  magnitudes,  were  found  and 
attributed  to  data  source  considerations. 

The  inter-annual  patterns  showed  a  cellular  structure  that  remains  in 
phase  through  all  the  analysis  changes.  The  EBBI  pattern  was  in  agree¬ 
ment  with  the  wave-like  equatorial  pattern  we  observed  in  the  divergence 
field.  Magnitude  variations  were  most  apparent  after  implementation  of 
GDAS  changes  in  May  of  1986.  The  patterns  in  the  EBBI  field  and  magni¬ 
tudes  over  the  8  year  period  were  very  stable  which  is  indicative  of  a  homo¬ 
geneous  data  set. 
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(lla)  Zonal  Mean  Divergence  1979-1986 


(lib)  Zonal  Mean  EBBI 


Figure  11.  Inter-annual  variability,  (a)  200  mb  zonal  mean  divergence,  con¬ 
tour  interval  3  x  10-6  s-1  Vertical  lines  are  times  of  GDAS  changes,  (b) 
EBBI,  contour  interval  6  ”  K. 


5.  Mean  Zonal  Circulations 

The  Walker  cell  is  a  system  of  zonal  overturning  motions  in  the  south 
Pacific  which  is  one  component  of  a  broader  system  of  zonal  overturning 
motions  in  the  tropical  belt.  Its  motions  ascend  in  the  western  monsoon 
region  and  descend  off  the  eastern  South  American  coast.  Bjerknes  (1969) 
used  a  small  tropical  data  set  to  show  the  existence  and  describe  the  struc¬ 
ture  of  this  circulation.  The  timeliness  of  Bjerknes’  work,  in  conjunction 
with  the  increased  availability  of  tropica!  data,  led  to  several  related  stu¬ 
dies. 

Krishnamurti  (1971a)  studied  the  tropical  east-west  circulations  using 
200  mb  velocity  potential  data  from  JJA  of  1967.  In  1973  he  collected  and 
analyzed  the  same  variables  from  DJF  of  1969  to  complete  his  description  of 
the  seasonal  variability  in  this  field.  He  observed  that  200  mb  east-west 
overturning  motions  came  from  the  same  monsoon  region  described  by 
Bjerknes  in  1969  and  described  it  as  a  southern  extension  of  the  main  east- 
west  flow  pattern  about  the  equator. 

Our  analysis  confirmed  Krishnamurti’s  observations  that  the  Asian 
monsoons  are  the  center  for  global  divergent  motions.  In  addition,  the  east- 
west  overturning  motions  over  the  Atlantic  were  strongly  influenced  by  con¬ 
vection  over  South  America  during  each  of  the  seasons.  The  isopleths  in 
JJA  (figure  5)  showed  comparable  magnitudes  between  Hadley  and  east- 
west  motions.  Krishnamurti  measured  the  intensity  of  east-west  circula¬ 
tions  by  calculating  (IE): 


where  y!  and  y2  are  the  southern  and  northern  limits  of  the  belt  and  IE 
varies  along  x.  We  chose  the  same  latitude  belts  used  by  Krishnamurti:  in 
DJF  yj  and  y2  are  15 0  S  and  15  °  N  and  in  JJA  they  are  0  *  N  and  30 0  N 
respectively. 

Figure  12  shows  the  mean  zonal  divergent  winds  for  the  two  major 
convective  latitude  belts  described  above.  Positive  (negative)  values  are 
eastward-  (westward-)  directed  200  mb  irrotational  velocities.  The  magni¬ 
tudes  appear  to  be  very  even  between  the  two  seasons  but  the  positions  of 
the  circulations  shift  to  the  west  by  about  60°  from  DJF  to  JJA.  The 
easterlies  are  associated  with  an  anti-clockwise  cell  in  the  vertical  plane  with 
a  rising  branch  over  the  monsoon  region  and  a  descending  branch  over 
Africa.  The  westerlies  over  the  Pacific  Ocean  are  associated  with  the  so 
called  Walker  circulation. 

These  observations  agreed  with  the  spatial  patterns  in  Krishnamurti’s 
figures  except  for  slightly  higher  magnitudes  in  our  patterns.  Because  his 
analysis  was  based  on  a  very  limited  number  of  observations  taken  over  one 
season,  it  is  highly  likely  the  discrepancies  are  found  in  the  data  sources. 

Figure  13a  is  the  monthly  mean  200  mb  zonal  variations  of  the  irrota¬ 
tional  velocity  field  averaged  over  8  years  and  the  latitude  belt  from  15  °  N 
to  15  °  S.  The  sign  convention  is  the  same  as  above.  We  have  highlighted 
lines  of  upper-tropospheric  convergence  (divergence)  found  at  the  junction  of 
Low-High  (High-Low)  patterns.  The  annual  variations  of  circulations 
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the  ’westerlies’  which  form  the  200  mb  component  of  the  Walker  cell  in  the 
south-central  Pacific.  The  200  mb  divergence  waivers  from  130  °E  to 
170 0  E  and  from  60 0  W  to  80 0  W.  The  central  line  of  200  mb  convergence 
waivers  from  80 'W  to  110°W.  The  most  significant  feature  are  the  varia¬ 
tions  in  the  western  Pacific  source  and  southeastern  Pacific  sink  regions. 


The  EBBI  and  Divergence  data  appear  to  be  in  agreement  on  the  place¬ 
ment  of  major  features,  such  as  lines  of  convection  and  upper-tropospheric 
divergence.  We  can  infer  a  major  center  of  convection  in  the  EBBI  field 
during  JJA  at  90 0  E  which  is  slightly  east  of  the  200  mb  divergence  line  in 
figure  13a.  Similar  displacements  in  other  divergence/ convergence  lines 
between  figures  13a  and  13b  can  also  be  seen. 

Figure  14a  displays  the  inter-annual  variations  of  the  mean  zonal 
divergence  for  the  same  latitude  belt  as  in  Figure  13.  The  inter-annual  pat¬ 
tern  remains  consistent  with  the  patterns  described  for  the  intra-annual 
cycle  with  one  small  anomaly  in  the  EBBI  field  at  150 ’E  in  late  1982  where 
a  minimum  in  the  band  of  positive  values  implies  strong  subsidence. 
Corresponding  to  this  minimum  is  a  belt  of  inferred  convection  from 
150  °W  to  120  °W  which  is  a  zone  which  normally  implies  subsidence. 
Closer  examination  of  the  divergence  field  shows  the  same  type  of  anomaly 
pattern  with  the  observation  of  a  maximum  indicating  a  reversal  from  the 
normal  easterly  flow  to  a  westerly  component.  Since  this  corresponds  to  the 
area  and  timing  of  the  1982-83  El  Nino  event  we  looked  at  this  period  in 
more  detail. 

Figure  15  is  similar  to  figure  12  except  it  shows  the  mean  zonal  flow 
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over  the  two  most  active  seasons.  The  January  1982  and  1983  data  are 
averaged  over  15  *  N  to  15  °S  and  the  July  1382  and  1983  data  are  averaged 
over  0  *  N  to  30  °  N.  The  figure  only  cover  the  longitudes  from  the  Indian 
Ocean  to  the  east  coast  of  South  America  to  highlight  the  Pacific  circula¬ 
tion.  There  is  a  slight  variation  in  magnitudes  between  1982  and  1983  but 
the  most  significant  contrast  is  the  60  °  and  30  °  iag  between  Januaries  and 
Julies  respectively.  These  lags  imply  movement  of  the  western  Pacific  source 
region  associated  with  the  ascending  branch  of  the  Walker  cell. 

In  Figure  16,  the  January  1982  convection  is  centered  over  the  PMC 
with  easterly  flowing  divergent  velocity  vectors  on  both  sides  of  the  equator 
with  convergence  of  Baja,  Mexico,  and  the  SPCZ  off  the  South  American 
coast.  In  January  1983,  the  convection  has  moved  to  165  °W  with  divergent 
velocity  maxima  decreasing  by  25%  from  1982  and  the  divergent  flow  in  the 
western  Pacific  has  reversed  to  flow  towards  the  west. 

In  Figure  17,  The  ITCZ  in  July  1982  is  clearly  marked  between  5 °  N 
and  10  °N  from  the  PMC  to  the  Baja,  Mexico.  The  western  ’source’  is  cen¬ 
tered  in  the  NH  at  150 0  E  with  a  maxima  in  the  EBBI  of  17.6  °  and  another 
maxima  of  21.4°  is  found  on  the  northwest  tip  of  Venezuela.  The  Walker 
cell  circulation  is  apparent  in  the  divergent  velocity  vector  field.  In  July 
1983  the  western  Pacific  convection  is  centered  at  120  °E  and  the  EBBI 
maxima  over  Venezuela  has  shifted  60°  northwest  of  its  1982  position.  The 
EBBI  infers  a  break  in  the  ITCZ  at  130  °E. 

In  summary,  the  analysis  of  this  data  generally  agrees  with  Ivrish- 
namurti  (1971a/b,1973)  who  also  found  stronger  velocities  in  the  JJA  than 
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DJF  cycle.  We  observed  the  seasonal  cross-equatorial  movement  in  the  mon¬ 
soon  convection  which  feeds  the  global  and  tropical  energy  systems,  and  the 
persistence  of  the  patterns  over  the  8  year  period.  We  also  observed 
significant  variations  between  January  1982  and  January  1983  and  between 
July  1982  and  July  1983  data.  The  1982  data  is  prior  to  the  start  of  the  El 
Nino  and  the  1983  data  is  during  the  warming  event. 
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Figure  14.  Inter-annual  :ycle  averaged  over  15  °  N  to  15  °  S  for  (a)  zona!  irro- 
tationa!  velocity  component,  contour  interval  3  m/s,  (b)  EBBI,  contour 
interval  10  a  K. 
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Figure  16.  Mean  EBBI  ( 3  K)  with  superposed  irrotational  velocity  vectors 
(m/s).  EBBI  contour  interval  8  0  K. 
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6.  200  mb  Divergence  and  EBBI  Correlations 

Kasahara  (1987b)  points  out  the  diabatic  NNMI  schemes  currently  used 
in  forecast  models  underestimate  the  release  of  latent  heat  during  the  first 
several  hours  of  the  model  run  and  this  leads  to  weakened  tropical  divergent 
motions.  Different  procedures  have  been  studied  to  improve  the  accuracy  of 
computing  tropical  heating  rates,  including  computations  from  satellite 
derived  data.  Richards  and  Arkin  (1981)  developed  an  algorithm  for  calcu¬ 
lating  tropical  rainfall  rates  from  infrared  flux  data  derived  from  satellite 
measurements.  They  were  able  to  correlate  regions  of  rainfall  with  low 
infrared  fluxes  and  latent  heat  release  can  then  be  computed  from  these 
rainfall  rates.  Kasahara  et.  al.  (1987a)  also  studied  this  problem  using 
fluxes  converted  to  EBBI  fields.  They  used  EBBI  as  the  difference  between 
260  °K  and  Te  because  it  is  the  approximate  value  of  equivalent  emission 
temperature  for  the  earth-atmosphere  system.  They  then  correlated  this 
field  with  the  u  field.  Their  results  showed  the  dividing  line  between  ascend¬ 
ing  and  descending  motion  is  approximately  258  ’  K. 

The  use  of  OLR  flux  as  proxy  data  for  conventional  observations  would 
partially  solve  the  data  gap  in  tropical  latitudes.  Convection  inferred  from 
OLR  flux  should  be  able  to  provide  a  gridded  divergence  field  in  these  lati¬ 
tudes.  In  order  to  achieve  this,  we  need  to  quantify  the  correlation  between 
the  two  fields. 

Julian  (1984)  and  Kasahara  et.  al.  (1987b)  are  two  examples  of  on¬ 
going  research  to  accurately  quantify  the  relationship  between  OLR  flux  and 
divergence.  They  have  had  some  success  based  on  very  limited  sets  of  daily 


observations.  In  the  present  study  we  use  spatial  and  temporal  correlations 
to  our  Tg  and  divergence  fields  to  see  if  significant  correlations  are  observed 
over  the  8  year  data  set.  The  methodology  used  in  computing  the  different 
correlations  can  be  found  in  the  Appendix. 

Figures  18,  19  and  20  are  spatial  correlations  (Type  I)  between  mean 
divergence  and  EBBI  fields.  The  mean  values  are  plotted  as  a  function  of 
longitude,  and  the  correlation  between  the  two  fields  is  indicated  in  the 
figure.  The  correlation  in  the  zone  from  5  °  N  to  the  equator  is  .84  which 
suggests  EBBI  explains  70%  of  the  variance  of  the  divergence  in  the  zone 
whereas  the  correlation  in  the  25  °  N  to  30  *  N  zone  is  only  .26.  The  overall 
correlation  between  the  tropical  belt  between  30  *  N  and  30  °  S  is  .77  explain¬ 
ing  about  60%  of  the  variance. 

Next  we  applied  temporal  correlations  (Type  II)  to  January  and  July 
fields.  The  results  are  shown  in  Figures  21  to  24  where  the  .40  and  .70  con¬ 
tours  of  the  correaltion  coefficient  are  shown.  The  .70  contour,  which  is  the 
significant  level  we  determined  for  this  set,  is  shaded.  Accompanying  the 
correlations  are  the  mean  EBBI  and  divergence  fields  for  comparison. 

Figures  21  and  22  show  similar  placements  in  the  inferred  divergence 
centers  over  South  America  and  the  PMC  and  the  same  elongated  figure  to 
the  southeast  of  the  PMC.  The  EBBI  pattern  infers  divergence  extending 
well  west  of  the  PMC  into  the  Indian  sub-continent  and  also  over  central 
Africa.  Neither  of  these  patterns  is  found  in  the  divergence  patterns.  There 
is  some  agreement  in  inferred  convergence  zones  centered  at  15  °N  160  °  W, 
over  the  Red  Sea  on  the  northeast  coast  of  Africa,  and  in  the  Caribbean 


v 


Sea. 
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The  results  on  correlations  are  not  as  positive  as  we  were  expecting.  At 
the  significance  level  computed  for  the  8  months  in  the  set,  there  are  only 
small  areas  of  significantly  tested  correlations  between  the  two  sets  of  data. 
In  the  small  areas  that  are  correlated,  the  most  significant  grouping  is  in  the 
vicinity  of  the  PMC  as  we  would  expect.  The  central  Pacific  convergence 
area  described  above  also  shows  some  good  correlations  as  do  other  conver¬ 
gence  zones  off  the  west  coast  of  South  America,  and  in  the  Caribbean  sea. 
It  is  surprising  that  in  these  nominal  areas,  convergence  areas  are  better 
represented  than  areas  of  divergence. 

Figures  23  and  24  are  similar  figures  for  July.  The  same  dominant 
divergence  centers  for  this  season  are  again  observed.  The  PMC,  the  Cen¬ 
tral  American  isthmus,  and  central  Africa  all  reflect  strong  200  mb  diver¬ 
gence.  There  is  a  nice  representation  of  the  inter-tropical  convergence  zone 
(ITCZ)  in  the  EBBI  field  that  is  not  as  well  represented  in  the  divergence 
field.  Only  small  correlation  patterns  appear  at  the  significant  level  and 
areas  of  convergence  are  again  as  well  represented  as  areas  of  divergence. 

Our  objective  in  this  section  was  to  quantify  the  relationship  between 
OLR  flux  and  the  200  mb  divergent  circulation.  We  did  this  by  computing 
several  correlations  between  our  independent  data  sets.  The  first  correlation 
of  band  average  mean  fields  in  Figures  18,19,  and  20,  supported  our  expec¬ 
tations  by  returning  high  correlations  in  belts  of  divergence,  weak  values  in 
belts  of  convergence,  and  overall  high  return  for  the  entire  tropical  belt.  The 
next  set  of  correlations  were  disappointing  because  only  isolated  patterns  of 
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significant  correlations  were  observed.  In  all  we  were  successful  in  showing 
latitude  belts  of  convection  are  well  correlated  with  OLR  fluxes  but  were 
less  than  successful  in  finding  significant  correlations  on  individual,  point- 
to-point  computations.  Because  these  values  are  needed  to  make  headway 
in  quantifying  the  relationship  between  OLR  flux  and  200  mb  divergence  we 
have  to  say  we  were  not  successful. 
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(22a)  Pattern  Correlation  for  January 
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(23a)  Pattern  Correlation  for  July 


7.  Summary 

We  studied  the  variability  of  the  200  mb  tropical  irrotational  wind 
field  and  its  relationship  to  variabilities  in  OLR  fluxes.  Horizontal  diver¬ 
gence  is  an  important  component  of  atmospheric  circulations  in  the  tropics. 
Until  recent  advancements  in  satellite  derived  data  the  measurement  of 
divergent  components  was  very  limited.  These  improvements,  along  with 
timely  studies  of  variations  in  tropical  circulation  by  Bjerknes  (1969)  and 
Krishnamurti  (1971a,  1971b,  and  1973),  were  important  to  focusing  research 
on  this  field. 

Velocity  potential  and  OLR  flux  data  for  the  8  years  from  January 
1979  to  December  was  provided  by  the  NMC  Climate  Analysis  Center.  The 
extent  of  this  data  allowed  us  to  examine  intra-  and  inter-annual  variabili¬ 
ties  over  an  extended  period.  (Bjerkne’s  and  Krishnamurti’s  work,  cited 
above,  were  based  on  very  limited  data  set)  Other  studies,  such  as  Kasahara 
et.  al.  (1987a)  and  Julian  (1984),  studied  the  problem  of  quantifying  the 
relationship  between  divergence  and  OLR  flux  fields.  We  also  examined  this 
problem. 

Our  first  observation  about  the  200  mb  tropical  divergent  velocity  field 
is  the  dominance  of  the  western  Pacific  source  region  associated  with  mon¬ 
soon  activity.  This  region  drives  both  the  meridional  and  zonal  overturning 
motions.  Smaller,  but  important  source  regions  are  found  over  South  and 
Central  America  as  well  as  Africa.  The  most  important  variation  we 
observed  concerning  these  source  regions  is  their  intra-seasonal  transitions 
to  the  summer  hemisphere.  A  second  observation  brought  out  by  seasonally 


averaged  data  was  the  variability  in  the  NH  flow  while  the  SH  had  only 
slight  variations.  This  hemispheric  asymmetry  is  in  response  to  differences 
in  zonally  asymmetric  forcing  features  between  the  hemispheres.  Our 
results,  in  general,  compared  favorably  to  those  of  Arkin  et.  al.  (1986)  and 
Lau  and  Boyle  (1987)  who  studied  similar  variations  in  tropical  circulations. 

Intra-  and  inter-annual  variations  in  the  zonal  mean  meridional  circula¬ 
tions  indicated  the  same  seasonal  meridional  transitions  found  in  the  sea¬ 
sonal  means.  In  addition,  the  ascending  and  descending  branches  were 
clearly  observed  in  the  monthly  mean  data.  Another  feature  was  the  strong 
annual  cycle  we  found  in  the  tropical  belt  from  30°  N  to  30°  S  and  the  sea¬ 
sonal  fluctuations  in  the  wind  patterns  driven  by  the  ascending  branch  of 
the  equatorial  Hadley  cell  are  quite  clear.  The  winds  shift  from  north  to 
south  in  6  month  periods  and  this  pattern  remains  spatially  and  temporally 
consistent  over  the  entire  8  years  with  only  magnitude  variations  associated 
with  changes  in  the  GDAS.  These  results  were  in  general  agreement  with  a 
study  of  zonal  mean  meridional  tropical  circulations  by  Oort  and  Rasmus- 
son  (1970). 

The  intra-  and  inter  annual  cycles  of  the  mean  zonal  EBBI  and  diver¬ 
gence  fields  showed  a  strong  relation  between  these  independent  data  sets. 
Seasonal  transitions  and  persistent  annual  cycles  are  reflected  in  both  sets. 
The  divergence  field  is  noisier  than  the  EBBI  field.  We  expected  this  because 
the  EBBI  is  very  homogeneous  in  comparison  to  the  divergence  data. 

Next  we  focused  on  the  divergent  motions  on  the  zonal  plane.  We 
found  strong  evidence  supporting  the  results  of  Krishnamurti  (1973)  that 
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the  monsoon  region  is  the  source  for  east-west  overturning  motions  in  the 
tropics.  In  particular  the  Walker  circulation  in  the  equatorial  Pacific  was 
strongly  supported  by  these  fields.  The  east-west  movements  of  the  mon¬ 
soon  caused  some  variations  in  the  spatial  positioning  of  these  motions. 
Again,  on  the  inter-annual  scale,  the  most  striking  feature  was  the  per¬ 
sistence  in  annual  cycles  from  year  to  year.  This  persistence  held  in  the 
qualitative  examination  of  the  divergence  and  EBBI  data. 

One  anomaly  in  the  zonal  mean  flow  which  proved  interesting  was  a 
noticeable  anomaly  in  the  inter-annual  patterns  of  the  zonal  and  EBBI 
fields.  This  happened  during  the  El  Nino  period  of  late  1982  and  1983  in 
the  central  equatorial  Pacific.  Further  investigation  brought  out  clear  varia¬ 
tions  between  patterns  prior  to,  and  in  the  middle  of,  this  event. 

Finally  we  computed  correlations  between  the  divergence  and  EBBI 
data.  Kasahara  et.  al.  (1987)  and  Julian  (1984)  both  have  studied  this  rela¬ 
tionship  with  only  limited  success.  We  were  able  to  find  a  strong  positive 
correlation  in  the  latitudes  close  to  the  equator  where  convection  is  very 
strong.  Conversely,  we  found  very  poor  correlations  as  we  moved  further 
away  from  the  equator.  More  specific  correlations  on  a  point-to-point  basis 
for  January  and  July  data  were  less  successful  because  only  small  areas  of 
significant  correlation  were  found. 

This  study,  based  on  a  relatively  larger  data  set  of  tropical  200  mb 
divergent  circulations  than  have  been  found  elsewhere,  confirmed  the  results 
of  earlier  works  based  on  limited  data.  This  study  also  contributed  to 
understanding  the  intra-  and  inter-annual  variations  in  these  circulations. 
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The  renewed  interest  in  these  circulations,  and  their  importance  to  general 
circulation  and  numerical  weather  prediction,  has  occurred  because  of  the 
wider  availability  of  data.  We  feel  continued  research  in  this  field  is  war¬ 
ranted  because  it  is  an  important,  yet  unsolved,  component  of  general  circu¬ 
lation  theory. 


Appendix:  Computational  Techniques 


1.  Central  Finite  Differencing 

The  irrotational  velocity  vectors  were  calculated  directly  from  the 
potential  field  using  central  finite  differencing. 


1  dx 

a  cos  <f>  d\ 

1  dx  _  1 

a  d4>  a 
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2.  Zonal  Averaging 

Zonally  averaged  fields  were  computed  as  follows: 

i  144 

VzG)  -  TH  £  vu  (j  -  1  •  •  •  73)  (15) 

i-i 

where  i  represents  the  144  longitudinal  grid  points  from  0  °  to  360 0 ,  and  j 
the  latitudes  from  90  °  N  to  90  "S,  each  with  2.5°  spacing  interval.  From 
this  definition,  the  zonal  average  Uz  equals  zero: 


U*  =  « - - - T  I  dX  -  0  (16) 

*  27T  a  cos0  J  d\ 

which  shows  no  east-west  momentum  can  accumulate  around  a  latitude  cir- 


3.  Divergence 
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The  equations  used  to  compute  divergence  V  2x  from  the  velocity 


potential  field  follow: 
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We  then  used  the  following  central  differencing  scheme: 
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Similar  calculations  can  be  made  for  -A  — -  to  get: 

o<p  d<p 


a  ax  cos0 

00  a^ 


=  [(x  cos 0)iij+l  +  (x  cos0),  J_J  -  2(x  cos0),  j  ]  (23) 


where  on  our  2.5  °  grid  has: 


(A<i)2  (AX)2 


4.  Correlation 

In  section  6  we  used  several  correlations  which  we  will  briefly  describe 
here  along  with  their  significance  test. 

a.  Type  I  used  to  compute  Figures  17,  18,  and  19. 

1.  In  these  figures  we  correlated  the  8  year  band-averaged  mean 
grids  for  the  divergence  and  EBBI  fields. 
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Let  Xjj  be  the  8  year  (96  month)  mean  over  the  entire  domain.  That  is: 

i  96 

Xj.  =  -  Sxiin  (24) 

96  ,jn  v 

i  =  l  ,...,25  (Latitude),  j  =  1,...,144  (longitude) 

The  tropical  grid  goes  from  30°  N  to  30  °  S  (1,...,25  at  2.5  °  spacing  interval) 
latitudes  and  from  0°  to  357.5°  (1,...,144  at  2.5°  spacing  interval)  meridi¬ 
ans.  Next  we  computed  the  band  average  for  the  appropriate  latitude  belt: 

AJ  =  */  "  £(*■  i)  (25) 

j  =  1...144 


Finally,  we  compute  the  correlation  between  the  Divergence  and  EBBI  fields 
using  the  general  correlation  equation: 
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where  A  and  B  are  mean  fields  computed  as  described  above,  and  A  and  B 


are  their  means: 


A  = 


1  144 
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144 
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2.  t  Test  for  significance: 


t  = 


rVN-2 

Vl  -  r2 


(28) 


where  N=144  and  we  get  t  at  the  95%  confidence  level  for  N-2  degrees  of 


freedom  from  tables.  Here,  1 95  U2  equals  1.66.  But  all  144  values  are  not 
independent  and  hence  it  is  not  obvious  what  the  degrees  of  freedom  for 


this  case  is  and  therefore  no  significance  test  was  computed. 


b.  Type  II,  used  to  compute  Figures  21  to  24. 

1.  In  these  figures  we  compute  temporal  correlations  based  on  the 

8  Januaries  and  the  8  Julies  from  1979  to  1986.  The  tropical  grid  used  in 

these  figures  was  reduced  to  the  latitudes  from  20 0  N  to  20  °  S  which 
reduced  the  latitude  index  to  i  ( 1  ,...,17)  First  we  computed  the  mean  Janu¬ 
ary  and  July  fields  then  the  8  anomaly  grids  by  subtracting  the  annual 
mean  from  the  individual  monthly  grids  and  finally  the  correlations: 

E  j)  (YiJn-yu)] 

r  ss  - 1 _  (09) 

"  VfEtx.i.-?,,)2  (E(Yin,-y;j)2 

2.  The  t  test  is  used  to  check  the  significance  at  the  95%  confidence 
level  with  N  =  8,  and  1 956  equal  to  2.447.  This  makes  .70  the  significant 
correlation  level. 
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